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Abstract 

\ A universe described by braneworlds is studied in a cyclic scenario. As expected such an 

■ oscillating universe will undergo turnarounds, whenever the phantom energy density reaches a 
\ critical value from either side. It is found that a universe described by RSII brane model will 

^ ^ ■ readily undergo oscillations if, either the brane tension, A or the bulk cosmological constant, A4 

pi I is negative. The DGP brane model does not readily undergo cyclic turnarounds. Hence for this 

• model a modified equation is proposed to incorporate the cyclic nature. It is found that there is 

I always a remanent mass of a black hole at the verge of a turnaround. Hence contrary to known 

I ■ results in literature, it is found that the destruction of black holes at the turnaround is completely 

I out of question. Finally to alleviate, if not solve, the problem posed by the black holes, it is argued 

I""!. that the remanent masses of the black holes do not act as a serious defect of the model because of 

Hawking evaporation. 

U '. 

1 Introduction 

CO ■ Cyclic universe has always been a burning topic in the field of theoretical cosmology, since it is expected to 

^ . avoid the initial singularity by providing an infinitely oscillating universe. However cyclic universe confront a 

serious problem of black holes(BHs). If the BHs formed during the expanding phase survives into the next cycle 
. they will grow even larger from one cycle to the next and act as a serious defect in an otherwise nearly uniform 

V/^ ■ universe. With the passage of time the BHs will occupy the entire horizon and then the cyclic models will break 

■ away. In this paper we investigate the possibility of an oscillating universe in two of the well known models of 

■ brane-world gravity, namely, RSII brane and DGP brane models. 

. Randall and Sundrum [T,'? proposed a bulk-brane model to explain the higher dimensional theory, popularly 

' known as RS II brane model. According to this model we live in a four dimensional world (called 3-brane, a 

^ ■ domain wall) which is embedded in a 5D space time (bulk) . All matter fields are confined in the brane whereas 

• j ■ gravity can only propagate in the bulk. The consistency of this brane model with the expanding universe has 

r> , given popularity to this model of late in the field of cosmology. 

■ A simple and effective model of brane-gravity is the Dvali-Gabadadze-Porrati (DGP) braneworld model 
[2 which models our 4-dimensional world as a FRW brane embedded in a 5-dimensional Minkowski bulk. It 
explains the origin of dark energy(DE) as the gravity on the brane leaking to the bulk at large scale. On 
the 4-dimensional brane the action of gravity is proportional to whereas in the bulk it is proportional 
to the corresponding quantity in 5-dimensions. The model is then characterized by a cross over length scale 

rc = ^'^'^^ ^^^^ gravity is 4-dimensional theory at scales a << rc where matter behaves as pressureless dust, 

5 

but gravity leaks out into the bulk at scales a >> Tc and matter approaches the behaviour of a cosmological 
constant. Moreover it has been shown that the standard Friedmann cosmology can be firmly embedded in DGP 
brane. 

In the context of BHs and phantom energy accretion on BH, it should be mentioned that Jamil et al [4] 
studied charged BHs in phantom cosmology. Jamil in j5] has shown the evolution of a Schwarzschild Black 
Hole in Phantom-like Chaplygin gas Cosmologies. Primordial BHs in phantom cosmology and accretion of 
phantom DE on BTZ BHs were also studied by Jamil et al in [51 [7] • Nayak in ^ investigated the effect of 
Vacuum Energy on the evolution of primordial BHs in Einstein Gravity. Paolis in [9] studied BHs in bulk 
viscous cosmology. In the context of cyclic cosmology, it should be mentioned that Saridakis in [10] studied 
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cyclic Universes from general collisionless Braneworld models. Cai et al in [TT] investigated cyclic extension 
of the non-singular cosmology in a model of non-relativistic gravity. Cai et al in [T2| investigated cyclic and 
singularity-free evolutions in a universe governed by Lagrange-multiplier modified gravity. Moreover Cai et al 
in [13] showed that gravity described by an arbitrary function of the torsion scalar, can provide a mechanism 
for realizing bouncing cosmologies, thereby avoiding the Big Bang singularity. Non-singular cyclic cosmology 
without phantom menace was also studied by Cai et al in ^14j . 

We intend to study the effects and consequences of phantom energy accretion onto BHs in a cyclic scenario 
of the universe described by DGP and RSII branes. Our motivation is to find out if there is any remanent 
mass of BH when it undergoes a turnaround in a cyclic scenario. Babichev et al PJ] has shown that BH 
mass decrease with phantom energy accretion on it. Hence the BH will disappear before the turnaround in 
an oscillating universe. But Sun 16 provided a mechanism which showed that in an universe described by 
modified Friedmann cosmology the destruction of BHs is totally out of question, as there is always a remanent 
mass of a BH facing a turnaround. In this paper our motivation is to testify the above fact for branc-world 
cosmology and find out the fate of a BH undergoing phantom energy accretion in an oscillating universe. 

The paper is organised as follows: In section 2 we discuss the mechanism of cyclic universe in RSII brane 
model. Section 3 deals with an identical mechanism for DGP brane model. In section 4, we present an argument 
regarding Hawking evaporation of remanent BHs. Finally the paper ends with some concluding remarks in 
section 5. 



2 Cyclic Universe in RSII Brane Model 

The novel feature of the RS models compared to previous higher-dimensional models is that the observable 3 
dimensions are protected from the large extra dimension (at low energies) by curvature rather than straightfor- 
ward compactification. In RS II model the effective equations of motion on the 3-brane embedded in 5D bulk 
having Za-symmetry are given by [HI [H H [H HOI 121] 

'•'^^Gm'^ = -A4gpi. -I- KlTfj_i, + KgH^^ - Ef,i, (1) 

where 

4^1 ^4 . (2) 

2 ' 6 

and 111 1 

n^l^ = T^aC + T"rf^l^ + g qf.l^Tal^lT"^ - ^ (4) 

and Efj_i, is the electric part of the 5D Weyl tensor. Here K5, A5, t^,^ and A4 are respectively the 5D gravitational 
coupling constant, 5D cosmological constant, the brane tension (vacuum energy), brane energy- momentum 
tensor and effective 4D cosmological constant. The explicit form of the above modified Einstein equations in 
flat universe for RSII brane are 

3H'^A, + 4p+!^ p^ + ^U (5) 
+ = A,-kIp~!^PP^^P^-^1( (6) 



A, = IkI(a, + ]:kI\'\ (3) 



and 



The dark radiation U obeys 



^4 



U + AHU (7) 
where p and p are the total energy density and pressure respectively. 
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2.1 Phantom Energy Accretion 

We consider an homogeneous and isotropic universe filled with DE fluid, with DE density p and pressure p. For 
an asymptotic observer the black hole mass, M changes at the rate of [T5] 

M = AttAM^ [p + p) (8) 

Here the overdot denotes the derivative with respect to cosmic time. Moreover it has been considered that 
G = c = 1. We consider an universe dominated by DE. The Friedmann equation for the expanding universe is 
given by 

- 'ip (9) 

Where H = ^ is the Hubble parameter. The conservation equation for DE is given by 

p + 3H{p + p)=0 (10) 

Phantom DE has the equation of state w — ^ < —1. We get that p cx a~^^^+™\ which shows that p increases 
with the expansion of the universe. 

In this section we will investigate the dynamics of cyclic universe in RSH brane model. The modified 
Friedmann equation for RSH brane model is given by equation ([5]). Putting U = 0, the equation becomes. 

In the expanding phase of the universe, the phantom energy density, p increases. From the above equation we 
see that at the turnaround the phantom energy density is given by. 



Pc 



-X±^/y-a (12) 



where a = 

From the above value we see that when A > 0, then A4 < 0, and when A < 0. then A4 > 0. A > This 
shows that either the brane tension or the bulk cosmological constant has to be negative, so that the universe 
undergoes a bounce, and a possibility for the cyclic scenario is evident. Here pc stands for critical density. After 
the turnaround the universe begins to contract. The reason behind this being, that in the contracting phase the 
non-phantom components of the universe increase and begins to dominate the evolution. In this phase, again 
when the dominant energy density reaches the critical value given by equation (|12p . a bounce occurs. Thus we 
get an oscillating scenario. We intend to study the variation of BH mass with phantom energy accretion around 
it, in this oscillating cosmological model. 



2.2 Scenario Before Turnaround 

Using equations dS]) and (fTU)) we get 



Before turnaround, we have 

Now substituting the above value of H in equation ([T3| we get 



^=^Wp(1 + 7Tt)+^ (14) 



dM D 



/p(l + ^) + A. 



-.dp (15) 



where D = and pc is given by equation (fT2)) . Now integrating equation (fT5|) we get 

l + i:>M,V2Alog 



M ^ , (16) 
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Here pi and Mi denotes respectively the phantom energy density and the black hole mass at the moment when 
the phantom energy density just begins to dominate the evolution of the universe. In general, pi <^ pc and 
Pi < p 1^ Pc- So using equation ([T2l) we obtain ^ — >■ 0. Hence from equation (|16[). we obtain 



, (17) 



1 + D^h^/2\\og 



At the turnaround, p — pc and hence the black hole mass at the turnaround becomes, 

Mc =^ ^ , (18) 

/ A+pc + ^«+2Ap,+pg 



l + LiM,V2Alog 



This shows that there is a remnant mass of the BH when the turnaround occurs. Hence this result is different 
from the result obtained by Zhang [35]. It is quite clear from the above equations that initially, through 
phantom energy accretion, the BH mass decreases, until it reaches the minimum value at the turnaround 
in the expanding phase. For Mj ^ Mp = G^2^ Mc becomes independent of Mi 

Mc ^ 7 ^ , (19) 

A+Pc + Va+2Ap,+pg 



DV2Alog 



2.3 Scenario After Turnaround 

After the turnaround as expected the the universe will contract and consequently the phantom energy density 
p starts decreasing. During the contraction it is obvious that the time derivative of the scale factor will become 
negative, and as a result H becomes negative. So, here we take the value of H as. 



Substituting the above value of H in equation ([T3|) we get, 

dM D 



Now integrating the above equation we get. 



M = ^ , (22) 

A+Pc + Va+2Ap^+pg 



1 + DMcV2Xlog 



where p < Pc- The above equation shows that as the universe contracts, the BH mass continue to decrease. 
When p Pc, then the BH mass is, 

Mf ^ , (23) 

A+Pc+-^a+2Ape+pg 



1 + DMcV2X log 



For Mi ^ Mp, using equation (fT9|) we find that the final mass of BHs is 

2i:>V2Alog 



Mf^^c^ ^ , (24) 

2 on./:TTi„„ M+P^+V"^+2APc+Pr 



Hence Mf is independent of Mi. 
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(Variation of BH mass in the expanding cycle(before turnaround Variation of BH mass in tlie contracting cycle( after turnaround" 




Fig. 1 Fig. 2 



Fig 1 : The mass of the black hole is plotted against the increasing density of phantom dark energy 
before turnaround. Other parameters are fixed at a = —5, X = 10, D = 1.5, Mi = 10000, pi = 0.0001 



Fig 2 : The mass of the black hole is plotted against the decreasing density of phantom dark energy 
after turnaround. Other parameters are fixed at a = —5, A = 10, D = 1.5, Mc = 100, pc = 100 



3 Cyclic Universe in DGP Brane Model 

While flat, homogeneous and isotropic brane is being considered, the Friedmann equation in DGP brane model 
is modified to the equation 

where H — ^ ia the Hubble parameter, p is the total cosmic fluid energy density and Tc = is the cross-over 

scale which determines the transition from 4D to 5D behaviour and e = ±1 (choosing = 8ttG = 1). For 
e = +1, we have standard DGP(+) model which is self accelerating model without any form of DE, and effective 
w is always non-phantom. However for e — —1, we have DGP(— ) model which does not self accelerate but 
requires DE on the brane. It experiences 5D gravitational modifications to its dynamics which effectively screen 
DE. 

Like the original Friedmann equation, and contrary to RSII brane equation, the DGP brane equation does 
not readily support an oscillating universe undergoing turnarounds. This fact is quite obvious from the equation 
(PSj) . We can see that there is no possibility of turnaround. Therefore we will propose a modified DGP brane 
equation, that can support a cyclic universe. Sun et al in |23| proposed a method, by which they were able to 
modify the Friedmann equation into a form that avoids the Big Rip singularity and gives a bouncing cosmological 
model. Their investigations led to the fact that the character of physics changes remarkably near the Planck 
scale. We know that a de Sitter universe with a cosmological constant, A is similar to a BH. Also it has a 
temperature, T ~ /J. Hence by conjecturing physics at the Planck scale, they actually modified the definition 
of Hawking temperature, and subsequently obtained the modified Friedmann equation. Sun in [TB] used this 
modified equation successfully to study the phantom energy accretion on BHs. Hence, taking a leaf out of their 
book, we proceed to modify the DGP brane equation as follows: The proposed modified DGP brane equation 




From the above equation, we see that the turnaround occurs for p — p'^. 
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3.1 Scenario Before Turnaround 

The expression for H for the expanding phase is given by, 



\ 3 4r2 2rc , 



1-^ 



Using equation ([T^ and (P7|) we get, 

M2 



£ + — + — 

3 ' 4r2 ^ 2r<: 



where = 



Now integrating the above equation we get. 



M Mi 



7( 



3pJ, arctan 



1-4 



arctan 



(27) 



(28) 



3A - p'^ + 2p, 



I 



( 



arctan 



V(Pc-pO(3^-P.), 



9^2 + 3B2 (p'^ _ 2p) + pp^ + 3A (p + - 3B 



V 



25 J3p^ (3A + p) (3i?2 _ 3A - p^) ( 1 - ^ 



- arctan 



9^2 + (p^ - 2p0 + p,p'^ + 3^ p, + p^ - 3^2 



V 



2B^3p^ (3A + p,) (352 _ 3^ _ p, ) _ 
/ 



^3A - 3^2 + p'^ 



arctan 



V 



1-4) \ 



^•iA - 3^2 + p'^ 



arctan 
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V 



^3A - 3^2 + p'^ 



(29) 



Where A = and B 



2ra 



Now as in the case of RSII brane, pi <C p^, which imphes that ^ — 0. Hence from the above equation we 



obtain, 



1^1 



D' 



3p^ arctan 



3A - p^ + 2p 
2v/(Pc-p)(3^-p). 



arctan 



3A-pi 



2\/3A^, 



/ 



35 Vp^ 



^352 - 3A - p^ 



arctan 



9^2 + 352 (p'^ _ 2p) + pp^ + 3.4 (p + p^ - 352 



V 



2BW3p^ (3A + p) (3B2 _ 3A - p^) 1 - f 



+- 



arctan 
arctan 



^2,A ~ 3B^ + p'^ 
At the turnaround p = p'^. Hence we get, 



V 



9/12 ^ 3^2^/^ _^ 3^ i^p,^ _ 3^2 
2B^9Ap'^ (3B2 _3A-p'J 



^3A - 352 + p^ 



arctan 
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^3A - 352 



(30) 



1 



1 



Mr. M,; 



/3p^ — arctan 



3A-P^ 
2^3^ 



6 



TT 

arctan 



: arctan 



v/3A - 3B2 + 



(31) 



The above equation gives the remanent mass of the BH at the turnaround. Hence we see that the mass of BH 
gradually decrease during the expanding phase of cyclic universe, due to phantom energy accretion and finally 
becomes minimum at the turnaround, when it is called the critical mass, given by equation ([30]). For Mi 3> Mp, 
Mc becomes independent of Af,; and is given by the following relation, 



Mr 



— y-ip'c ~ arctan 



3A-P'c 



^3ij2 -3A-P 

6ByM 



arctan 

2 



^3A - 3^2 + p'^ 



arctan 



2V3m; 

9^2 + 3B^p'^ + 3A {p'^ - 3^2) 
2B^9Ap'r (3S2 -3A- p'J 

/7c 



y/3A - 3B2 + p'^ 



(32) 



3.2 Scenario After Turnaround 

Just like the RSII model, in DGP model as well, the universe starts to contract after turnaround, and the 
non-phantom components starts to dominate the evolution. The expression for H is given by. 



H 



1 



3 4r2 



Using equation and ([551) S^t, 



dM 



1 



D' 



1 

4rJ 



Integrating the above equation we get. 



M Mi 



0' fa 



'3p'^ I arctan 



3BJp'^ 



^3^2 -3A-p'^ 
I 



3A-p'r + 2p 
V(Pc-p) {iA-p)^ 

I 



1-4 



— arctan 



3 A - p'r + 2p, 



arctan 



V 



V(Pc-pO(3A-p,:). 



9^2 + 3^2 {p'^ - 2p) + pp'^ + 3A{p + p'^- 3B 



V 



2BJ3p', {3A + P) {3B^~3A-p',)[l 



I 



■ arctan 



9^2 + 3^2 (p^ - 2pi) + p.p'r + 3^ 



3^2 



V 



2BW3p^ {3A + pO (352 -3A- p'^ 



1-4 



^3A - 3B2 + p'^ 



arctan 



V 



^3 A - 3^2 + p'^ 



arctan 



7 



^3A - 3B2 + p'^ 
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(33) 



(34) 



(35) 



Now as in the previous section, p^ ^ p^, which implies that ^ — ^ 0- Hence from the above equation we obtain, 
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M Af, 



farctan ( ^^^-P'^ + ^P 



arctan 




2B^ip', (3A + p) (352 _ 3^ _ _ ^) 

9^2 + 3S2p^ + 3A (p^ - 3^2^ 
2B^QAp'^ (3B2 _3A-p'J 



+ 



V'3A - 352 + 



V 



( ( M-^ ] 

arctan 



y/3A - 352 + 



arctan 
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At the turnaround p = p'^. Hence we get, 



Mr Af, 



arctan 



3A-p'r 



35, 



V352 - 3A - y2 



arctan 



9^2 + 352p^ + 3A (p^ - 352 
25V9M(352- 3A-p', 

7c 




(36) 



, I arctan ( — | (37) 

^3A - 352 + p^ yv^3A-352 + p^yJ 

As before the above equation gives the remanent mass of the BH at the turnaround. For Af, ^ Afp, Mc becomes 
independent of Mi and is given by the following relation, 



arctan 



3^-p'c 



35, 



v/352 - 3A - p^ \ 2 



TT 

arctan 



65, 



v/3A - 352 



; arctan 



2^33^; 

9^2 + 352p^ + 3A (p'r 



2B^%Ap'^{iB^ -iA- p'^ 




^ZA - 352 



(38) 



4 Hawking evaporation of black holes 

From the above evaluation it is clear that there is always a remanent mass of the BHs at the turnaround in an 
oscillating universe. So contrary to existing literature we conclude that there is no possibility of destruction of 
BH . The BHs formed during the expanding cycle of the cyclic universe survive into the next cycle and eventually 
grow in size. Therefore they create undesired non-uniformity in a nearly uniform universe. Eventually the BHs 
will occupy the entire volume of the horizon and will be responsible for the destruction of the cyclic models. This 
is a serious defect indeed! Hence the problem posed by the BHs in a cyclic universe still stands un-eliminated. 

But in [24 it has been argued that for a BH with mass M — 10^ Afp, Hawking evaporation takes place in 
time r • 



25ttM^ 
Mi 



10 sec and ultimately the BH becomes non-existent, thus causing no problems. Here Afp 



represents Planck mass. In our above calculations we have considered G 
for RSH brane we get, 

Af ^ ~ 



Afp 2 = 1. Rewriting equation 



DMiV2X\og 



TIT 



(39) 
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D is taken as a constant of the order unity and log ( '^p^^v'^^ — p^+Pc \ ^ ^2 ^j^gjj j^ave Mc ~ . 

For a BH with Mi ^ Mp, Afc ~ Mp, i.e., the remanent mass of the BH is of the order of Planck mass. 
Hence the remanent BH undergoes Hawking evaporation in time r ~ 10"^'^ sec at the order of Planck time. 
A similar evaluation is possible for DGP brane model as well. So fortunately, the remanent BHs do not cause 
any problems. This gives a possible solution of the BH problem in cyclic universe described by brane-world 
scenario. 



5 Conclusion and Discussion 

A serious problem is posed by the existence of black holes in an oscillating universe. In [SI] Brown et al 
suggested that in an oscillating cosmology the black holes keep losing mass due to phantom energy accretion 
before totally disappearing before the turnaround. Babichev et al in [15] devised a successful mechanism which 
was in accordance with the result given by Brown et al. In this paper we have investigated the outcome of 
phantom energy accretion on black holes in a cyclic universe described by brane-worlds. It is seen that RSII 
brane model readily incorporates the oscillating nature in its framework. The only condition being the negativity 
of the brane tension, A or the bulk cosmological constant, A4. But unlike RSII brane model, DGP brane does 
not readily support the cyclic nature of the universe. So a modified DGP brane equation, that supports cyclic 
turnarounds has been proposed for our evaluation. 

It is found that for RSII brane model, during the expanding phase the black hole mass gradually decrease 
with the increase of phantom energy density and finally reaches a critical value at which the turnaround occurs. 
This result is absolutely consistent with the known results in literature. In the contracting phase the black 
hole mass again decreases with the increase in the non-phantom components of the universe. This is however 
contrary to our expectations. So it is understood that our evaluations after the turnaround are not really 
rigid. We see that the black holes in a cyclic universe reaches a remanent mass, Mc before turnaround. So, 
the remanent mass implies that the destruction of black holes is not a real possibility in the cyclic cosmology, 
with phantom energy turnarounds, for a universe characterized by brane gravity. However fortunately we find 
that the remanent masses of black holes at turnaround do not cause problems. The reason being that these 
remanent black holes Hawking evaporate in a time r ^ 10"^"^. 
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